Distribution of Velocity in a Viscous Fluid, 205 

the value 0'46 is more accurate, this being the average of the values obtained 
in the most satisfactory experiments. 

The experiments were all carried out in the Engineering Laboratory of 
Cambridge University, and the author is indebted to Professor Ewing, F.RS., 
for much kindly advice and encouragement given during the course of the 
research. 
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1. Introduction. — An expression for the velocity at which stream-line 
motion breaks down in cylindrical pipes has been obtained by Osborne 
Reynolds,* and together with many others he has measured the fall of 
pressure occurring in different lengths. Lord Kelvin,t Lord Rayleigh,| and 
Reynolds§ have investigated the stability of different kinds of flow for 
viscous and non-viscous fluids, but how nearly the theoretical conditions of 
velocity distribution actually occur has not hitherto formed the subject of 
research. Experiments have often been made on the variation of mean 
linear velocities in the eddying state, but when the motion is irrotational, 
these are too low to admit of accurate measurement by any method so far 
employed. 

The objects of the present paper are to determine these velocities, to 
investigate the change which takes place at the critical velocity, and to find 
the relation between the velocity and the pressure to which it gives rise in 
a " Pitot " gauge of the form used. 

The research has become possible chiefly through the aid of an exceedingly 
delicate pressure gauge recently designed by my colleagues Professor 
A. P. Chattock and Mr. J. D. Fry. I am indebted to both these gentlemen 
for much valuable advice. 

^ ' Phil. Trans.,' 1883, Part III, p. 935. 

t ' Phil. Mag.,' August, 1887. 

X * Phil. Mag.,' July, 1892, etc. 

§ ' Phil. Trans.,' A, 1895, vol. 186, p. 123. 
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Osborne Eeynolds has shown that in parallel channels there are two 
limits to stream-line flow. The lower one is the velocity at which motion, 
originally eddying, will, when left to itself, settle down into stream-line 
motion, and the upper limit is that at which the motion of a fluid, originally 
at rest, traversing a pipe m an unstable stream-line condition, suddenly 
breaks down into the eddying or turbulent state. 

The variation of the lower limit with change of temperature has been 
further investigated by Coker and Clement,* and similarly that of the upper 
limit by Barnes and Coker.f The present paper deals with the lower limit 
of stream-line flow, and the term " critical velocity " is used in that sense 
throughout. 

2. Description of the Apparatus, — In an earlier apparatus the pipe was 
about 2*5 cm. in diameter and the Pitot tube 20 cm. long. There was 
danger of the Pitot tube bending, the velocity distribution was not 
symmetrical, and the temperature effects were considerable. 

In the new apparatus the pipe was of glass, 5*088 cm. in diameter. It 
was arranged with a bell-mouth in the entrance chamber, so that the water 
should enter under circumstances conducing to steadiness of flow, and by 
means of a fine tube a stream of colour could be admitted to indicate the 
nature of the motion. The water was supplied from the hydraulic installa- 
tion at University College, Bristol, and there was no other discharge from 
the mains while the experiments were in progress. 

The Pitot and static head tubes are shown in fig. 1, D being a brass pipe 
in continuation of A. To insure accuracy in the motion of the Pitot tube, 
p, it was carried by a rigid rectangular framework composed of a tube ^^, a 
rod rr, and two cross-pieces E and F. The guides in which this frame slides 
are two short tubes G and H and the longer one KL, and to these are 
attached the cross-piece M and scale S. 

Thus, by turning the screw E, the Pitot may be traversed across the pipe, 
its exact position being indicated by the reading on the scale opposite an 
index on F. 

The static head tube is at Q, a hole being drilled in the side of the glass 
pipe. A connection to qq would, therefore, give the pressure due to the 
velocity and static head combined, whilst one to Q would give the static 
head only. 

The two passages from the Pitot tube pass to one end of a reversing 
device, the other end of which is in communication with the static head tube. 
The two sides of this circuit are connected to the limbs of the pressure 

^ ' Phil. Trans.,' A, vol. 201, pp. 45—61. 
t * Koy. Soc. Proc.,' v^ol. 74, p. 341. 
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gauge. By this means (1) the static head can be brought to one side of the 
gauge and the Pitot tube to the other ; (2) the arrangement may be reversed ; 
or (3) the gauge may be short circuited on itself. 

It was found useful to be able to damp the slight oscillations occurring 
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at some velocities in the pressure gauge. For this purpose a small ait 
chamber and two constricted passages were provided. They could be cut 
out of the circuit, and were only used for verifying previous observations. 
Mg. 2 is a general view. It will be seen that the discharge took place 
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into a tank at the end of the pipe, and from there the water passed to waste 
by an overflow. 

The whole apparatus is carried on a stout timber base, the pipe and its 
attachments at each end being very carefully centred so as to ensure a 
maximum and symmetrical flow. 
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3. Measurement of pressure, discharge, and tem;perature. — The pressure 
gauge is of the Chattock-Fry differential type, and in construction is very 
similar to that in use at the National Physical Laboratory.* As employed in 
fchis research, however, it depended for its action on the difference in density 
of two liquids. 

From fig. 3 it will be seen to consist essentially of a U-tube, of which Ih 
are the arms and aa the horizontal connection. The two parts of aa meet 
in the chamber 0, and the fluid pressure is applied at CO. The tubes aa 
and the lower portions of Ih contain carbon-tetra-chloride, whilst the upper 
portions of hb and the chamber are filled with water. A meniscus 
surface of separation is shown at s, and the manipulation of the gauge con- 
sists in keeping the surface, viewed through the microscope M, absolutely 
unaltered in position. This is done by turning the graduated disc of the 
screw, and so tilting the gauge through the small angle necessitated by the 
variation of pressure difference. There is in this way no motion of the 
liquids relatively to the glass tubes which contain them, and hence errors 
due to surface-tension, refraction, or viscosity are entirely avoided. 

The pitch of the screw was 1 mm., and on the circumference of the disc 
were 200 divisions. These could be further divided by a vernier, but its use 
was abandoned as being within the limit of sensitiveness of the meniscus 
surface and miscroscope. 



* See Stanton, * Proc. Inst. C.E.,' vol. 156. 
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A and B being the dimensions shown in fig. 3, and E the reading of the 
screw, the pressure difference is given by 

A 

JL Xvir« Of 




A/B was found to be 1*038 and (p'— p) was taken from measurement on a 
Westphalia Balance to be 0*600 at 17° 0. 

To obtain the mean velocity of the water, it was necessary to measure the 
discharge. This was done by collecting the water passing through the 
apparatus in a given time, and weighing it on a carefully-calibrated spring 
balance to the nearest 10 grammes. The temperature was obtained by a 
thermometer placed in the end tank close to the outlet of the pipe. 

4. MetJiod of making the Experiments, — In determining the velocity distri- 
bution, experiments were made with the Bitot tube in various positions, 
commencing near one side of the pipe and continuing at successive points 
across the diameter, the water running continuously while a complete series 
of observations was made. The temperature was noted, and the discharge 
taken as the mean of six independent measurements. 

The following was the general system adopted for reading the pressures : — 
For any position of the Bitot tube, connection with the pressure gauge was 
made by opening the taps, the gauge was adjusted, and the reading taken. 
The connections were then reversed, and the reading taken. To eliminate 
errors due to friction or surface-tension, the observation was then repeated by 
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Fig. 4. 



closing tlie tap S, setting the pressure gauge to the previously obtained 
reading, reversing the connections, and then re-opening S. This was con- 
tinued until there was no apparent motion of the meniscus surface. The 
difference between the two readings was then double that due to the velocity 
of the stream at the point. A check on the accuracy of the gauge was applied 
during each series of readings. The instrument being adjusted when short- 
circuited, its zero reading was taken. The pressure due to the velocity was 
then read on each side of this zero. 

In the earlier experiments errors were discovered in this way, and found 

to be due to a very slight eccentricity 
of the point of the screw. Tests were 
also made for displacement of the 
surface when no water was flowing, 
and occasionally a very small correc- 
tion had to be allowed. This was 
attributed to slight differences of tem- 
perature in the passages leading to 
the gauge. 

5. Determination of the Pilot Con- 
stant, — If a tube be placed in a stream, 
as shown in fig. 4, and it be assumed, 
as in the Bernoulli Theorem, that the flow is steady, and that the velocities 
of consecutive stream-lines are not widely different, we may treat the fluid 
immediately at the mouth of the tube as at rest, and apply the theorem as 
though unimpaired by viscosity. 

The pressure ;p in the stationary fluid is then given by 

V = ip'^^l 
for this is the pressure in the adjoining stream-tube, where the velocity is 
indefinitely small. The formula requires experimental verification. Since 
the pressure must be proportional to pv^ we may in all cases assume 
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If h is the height above that due to the static head, to which the fluid will 
rise, the velocity of the stream approaching is 

where h is the Pitot constant to be determined. This may be done by a 
comparison of the gauge readings with the discharge from the pipe. The 
total flux Q is 

27rj rvdr, 

Jo 
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where v is the velocity distant r from the centre, and a the radius of the pipe. 
And, since v=:e^ (E) where E is the gauge reading and c a constant, we may 
write 

fa. 


To determine the integral, curves were plotted having r as abscissae and 
r^(R) as ordinates. The areas enclosed by the curves were then calculated 
by the method of ordinates, and the values obtained used for finding c, and 
hence the Pitot constant k An example is given in Table I, in which 
readings were taken at 19 different points in the diameter of the pipe. 











Table I. 












r. 


2-25. 


2 -00. 


1-75. 


1-50. 


1-25. 


1-00. 


0-75, 


0-50. 


0-25. 


0-00. 


V(R) X 102 
r^(E)xl02 


4-24 
9-54 


7-07 
14 -84 


9-53 
16 -68 


11-21 
16-81 


12-08 
15-10 


12-72 
12-72 


12-99 
9-74 


13-26 
6-63 


13-44 
3*36 


13-52 
0-00 



Q = 59-8 CO. per second. ^ = 1 -05. 

The second line in the table contains the square roots of the means of 
corresponding readings on each side of the centre. In fig. 5 the r-y/(E) 
curve shows the relative flux through different elemental tubes, into which 
we may suppose the stream to be divided. Its area was obtained from the 
mean of 83 ordinates. The ^/(R) curve is proportional to that of velocity 
distribution. The mean value of Jc obtained from the six experiments 
described in the next section is 1*04, and in calculating velocities I have 
adopted this constant throughout. The number may be compared with those 
recently found for air, namely, 1*03, by Dr. T. E. Stanton* and 0*974 by 
Mr. E. Threlfall, F.E.S.t 

6. Velocity Curves for Different Bates of Discharge, — Experiments were made 
at six rates of discharge, the curve of distribution of velocity being obtained 
in each case, and the result checked by comparing the mean velocity from the 
Pitot , readings with that obtained from the actual measurement of the 
discharge. 

The critical velocity, as determined by the point at which the law of 

velocity distribution changed, occurred when the mean was about 4 cm. 

per second (this being, of course, the velocity at which stream-line motion 

ceased to be the stable flow), and above this there was a decided change 

in the way in which the velocities increased with an increasing discharge. 

* * Proc. Inst. C.E.,' vol. 156. 
t * Proc. Inst. Mech. Eng.,' 1904. 
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The observed readings were reduced to velocities as already described, and 
these are given in Tables II and III. The r^y(B) curves were plotted, and 
the values of k determined. 



Fia. 5. 
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The temperature varied between 17^ and 19° C. 



Table II. — ^Experiments below Critical Velocity. 



Mean 
velocity in 
centimetres 
per second. 




Velocity in centimetres per 


second. 




Value of 




Radii in centimetres 

^ 


• 




0. 


0-5. 


1 00. 


1-50. 


2-00. 


2-25. 


1-4:35 


2-43 


2-39 


2-13 


1-85 


1-33 


0-78 


1-05 


2-196 


3-56 


3-45 


3-22 


2-88 


2-16 


1-17 


1-03 


2-955 


4-77 


4-67 


4-48 


3-96 


2-49 


1-50 


1-05 
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Table III. — Experiments above Critical Velocity. 



Mean 
velocity in 
centimetres 
per second. 


Velocity in centimetres per second. 


Value of 

k 


Eadii in centimetres. 

A . 


0. 


0-5. 


1-00. 


1-50. 


1-75. 


2-00. 


2-25. 


4-100 
5-373 

8-075 


6-43 

7-67 
10-57 


6-39 

7-51 

10-42 


6-12 

7-44 

10-25 


5-50 
6-93 
9-94 


5-14 
6-42 

9-38 


4-07 
5-53 

8-42 


2-48 
3-70 
5-30 

- 


1-04 
1-02 
1-05 



The curves in %. 6 are obtained by plotting the numbers in the last two 
tables. They show, as ordinates, the actual velocity of the stream at 



Fig. 6. 
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different distances from the centre of the pipe, the figures on the curves 
denoting the mean velocity as calculated from the discharge. 

It will be seen that below the critical velocity, the curve approximates to a 
parabola, and that afterwards there is a gradual change in form in the direction 
of a uniform velocity in the greater part of the pipe with a sudden fall 
to zero at the walls. These results are contrary to the frequent statement 
that, when the conditions of flow change from stream-line to eddying 
or rotational motion there is an equally sudden change in the law of distribu- 
tion of linear velocity from the parabolic to nearly uniform velocity over the 
cross-section. 

Fig. 7. 
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There appears to be no slipping at the walls of the tube. The nature of 
the change occurring at the critical velocity is exhibited in fig. 7, in which 
abscissae are mean and ordinates corresponding actual velocities at the 
different radii marked on the curves. From these it appears that, below the 
critical, the velocity at any fixed point is proportional to the total flux, and 
may be expressed approximately by the formula 

r being the distance from the axis of the pipe, v^, the mean velocity, A the 
ratio of maximum to mean, and B and n constants. In this case 

A = 1-62, B = 0-0815, n = 3-2. 
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According to the mathematical theory of viscosity* the surface graphically 
representing the distribution of velocity is a paraboloid, and the equation is 
of the above form, but having the constants 

A = 2, B = 2la\ n = 2. 

The two cases are, however, not strictly comparable for, whilst the theory 
takes no account of terminal conditions, in the experiments these are of 
importance. There would probably be accelerations and retardations over the 
cross-section, tending to equalise the velocities and the tubes of flow would 
not then be strictly parallel to the axis. It is, I think, evident that the 
variation of the curves from the parabola is in part due to these causes, and 
in part due to the fact that any resistance or obstruction across the pipe also 
tends to equalise the velocities. Beyond the critical velocity and within the 
very narrow limits of these experiments the curves in fig. 7 are approximately 
parallel when not near the sides of the pipe. At these low velocities, too, the 
distribution is represented roughly by the equation 

,8 



V =: 2'20-0'36r^"^-h {V04c-~^]vq. 



1750. 



In practical calculations in Hydraulics it is often assumed that the tangential 
resistance per unit area in a pipe is 

where / is a numerical constant depending on the nature of the surface. 
Accepting this as an expression of observed facts, Professor Lamb deduces an 
interesting conclusion.f If v denote the mean velocity at any point, we have 
at the surface 

u dvldr = \fpVf^, 

N"ow let I be the radial distance between two surfaces moving with relative 
velocity Vq in the regular " laminar " flow which would give the same 
tangential stress, and we find 

v^l = 2^1 1 fp. 

In the experiments, the last is the only one to which this formula can be 
applied. Thus taking fjujp = 0*0103, / = 0*005, and v^ = 8*075, we have 
I = 0*51 cm. The point at which the velocity is equal to the mean is found 
from fig. 6 to be almost exactly 0*5 cm., an agreement which is remarkably 
close for an approximate method. 

Summary and Conclusion, — The results may be briefly summarised as 
follows : — 

■^ Lamb's * Hydrodynamics,' p. 521. 
t Ibid., p. 574. 

Q 2 
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1. The experiments provide a partial confirmation of the theoretically 
obtained law of velocity distribution, but show that this distribution can only 
be obtained under very special conditions, of which absolute freedom from 
obstructions and end effects are important ; and hence 

2. When the flow is direct and stream-lines exist, the velocity distribution 
is not necessarily exactly that which may be described as characteristic of 
" normal " flow. 

3. At the critical velocity the irrotational straight line motion ceases and 
is followed by one in which the paths of the particles of fluid are eddying and 
turbulent. The law of distribution of mean linear velocity parallel to the axis 
simultaneously changes from the parabolic (or approximately parabolic) to that 
typical of eddying motion. 

4. The critical velocity in question (being that at which eddying motion 
ceases to be transformed into direct motion, and not that at which a highly 
unstable stream-line motion is suddenly disturbed), is not accompanied by a 
sudden change in the velocity parallel to the axis at any point in the cross- 
section. On the other hand, as the total flux increases, the experiments show 
a gradual transition from one state to the other, due to the change which has 
occurred in the law of velocity distribution. 

5. The observations have little bearing on the upper limit of stream-line 
flow, as observed by colour bands. They indicate, however, that the unstable 
direct motion would follow an approximately parabolic law of velocity 
distribution (as represented by the equation obtained for stream-line motion) 
and that at the higher critical velocity this distribution would suddenly 
change to that represented by the equation given for eddying motion. In this 
case then, instead of a gradual change of velocity, there would actually 
be sudden and large changes in the velocity parallel to the axis at different 
points in the cross-section of the pipe. 

6. The " Pitot Law " {v = ^2gh\ is at least approximately true at 
exceedingly low velocities. 

In conclusion I must mention my indebtedness to Professor E. M. Ferrier 
for his uniform kindness during the progress of the work, and for allowing 
portions of the apparatus to be constructed in the workshops of the College. 



